Extraordinary magnetoresistance of organic semiconductors : Hopping conductance 

via non-zero angular momentum orbitals 
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Highly-anisotropic in-plane magneto-resistance (MR) in graphite (HOPG) samples has been re- 
cently observed (Y. Kopelevich et al, arXiv:1202.5642 ) which is negative and linear in low fields 
in some current direction while it is giant, super-linear and positive in the perpendicular direction. 
In the framework of the hopping conductance theory via non-zero angular momentum orbitals we 
link extraordinary MRs in graphite and in organic insulators (OMAR) observed in about the same 
magnetic fields. The theory predicts quadratic negative MR (NMR) when there is a time-reversal 
symmetry (TRS), and linear NMR if TRS is broken. We argue that the observed linear NMR could 
be a unique signature of the broken TRS both in graphite and organic compounds. While some local 
paramagnetic centers are responsible for the broken TRS in organic insulators, a large diamagnetism 
of our HOPG samples may involve a more intriguing scenario of TRS breaking. 
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I. INTRODUCTION 

Transverse resistance, R(B), of isotropic conductors 
with Bloch electrons increases with an applied magnetic 
filed so that MR = R(B)/R(0) - 1 is positive and 
quadratic in B 0). A large positive MR is also used to 
be a hallmark of the hopping conductance in doped insu- 
lators [4j. While the carrier s-wave function on a donor 
(or acceptor) is spherically symmetric in the absence of 
the magnetic field, it becomes cigar-shaped squeezed in 
the transverse direction to the field [5-9]. This leads to 
a significant decrease in the overlap of the wave-function 
tails of two neighboring donors, and hence to a signifi- 
cant increase of resistivity (positive MR), which is also 
quadratic in low magnetic fields. 

However there are exceptions which do not show these 
canonical MRs. For instance some semiconductors and 
semimetals (e. g. bismuth), where open Fermi surfaces 
are unfeasible, show positive but linear MR. One of the 
theoretical possibilities for such a phenomenon is a so- 
called quantum magnetoresistance in semimetals having 
pockets of the Fermi surface with a small or even zero ef- 
fective mass (as the Dirac fermions in bithmuth, graphite 
and graphene), which might be in the ultra-quantum 
limit at rather low magnetic fields |10| . 

Also there is anomalous negative MR (NMR) observed 
in some hopping systems, for instance in amorphous ger- 
manium and silicon. Originally it has been attributed 
to magnetic- field dependence of spin-flip transitions be- 
tween sites when some fraction of them has a frozen spin 
plj . and/or to an increase of the density of localised 
states due to the Zeeman energy shift, \ibB [12 ]. This 
NMR is used to be small (much less than 1%) even in 
relatively high magnetic fields of about 1 Tesla. There 
are other theoretical mechanisms of NMR, in particu- 
lar the weak localization gives NMR which is often al- 
most linear in a certain field range. Such NMR smoothly 



evolves from a sub-linear magnetic field dependence at 
lower temperatures to super-linear field dependence at 
higher temperatures. A strong NMR exists in the clas- 
sical two-dimensional electron gas due to freely circling 
electrons, which are not taken into account by the Boltz- 
mann app roach. It is parabolic rather than linear at low 
fields [13( . The parabolic orbital NMR has been also pre- 
dicted by the gauge theory in two-dimensional strongly- 
correlated doped Mott insulators p4| . 

More recently a linear NMR has been observed in the 
longitudinal c-axis inter-layer current in the normal state 
of cuprate superconductors and in graphite in high mag- 
netic fields assigned to bipolarons in the former case [l5| 
and to a growing population of the zero-energy Landau 
level of quasi-two-dimensional Dirac fermions with the 
increasing magnetic field in the latter case Also 
a giant transverse NMR of nearly 100% was observed 
at low temperatures, with over 50% remaining at room 
temperature in graphene nanoribbons and attributed to 
some derealization effect under the perpendicular mag- 
netic field [17] . 

Puzzling magnetoresistance effects in several different 
7r-conjugated polymer and small molecular thin film de- 
vices have been observed named organic magnetoresis- 
tance (OMAR) pi]. OMAR reaches 10% at fields on 
the order of only 10 mT, and can be either positive or 
negative, depending on operating conditions. The effect 
is independent of the sign and direction of the magnetic 
field. The Zeeman energy does not account for the ob- 
served OMAR at ambient temperatures since it is too 
small, hbB ps 10 mK, in the field of 10 mT. The obser- 
vation of OMAR in hole-only devices [l8|, [l9| indicates, 
on the other hand, that the effect is hardly compatible 
with an exciton-based mechanism. 

An alternative model involving spin-dependent bipo- 
laron formation in deep potential wells has been proposed 
as a possibility 119]. In this model bipolarons (i.e doubly 
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occupied sites) block single polaron transport through 
bipolaron states causing positive MR, while an increase 
in polaron population at the expense of bipolarons with 
increasing magnetic field might cause negative MR, if 
the long-range Coulomb repulsion around each carrier is 
sufficiently strong. However, one might expect that this 
mechanism should depend strongly on the carrier density, 
whereas OMAR is only weakly dependent on current den- 
sity [Hj]. Finding a convincing explanation of OMAR is 
crucial for a better understanding of charge transport in 
organic semiconductors, actively used in light-emitting 
diodes, photovoltaic cells, field-effect transistors, and in 
spintronics [2ll |22| . 

Recentl y a n unusual orbital MR in HOPG has been 
observed [23(. In some current direction it is negative 
and linear in low fields with the crossover to the positive 
MR at higher fields, while in the perpendicular current 
direction MR is giant, super-linear and positive. In this 
paper we compare OMAR and graphite MR and pro- 
pose an explanation of both MRs in the framework of 
the hopping magneto-conductance via non-zero angular 
momentum orbitals 1241]. 



gamma-function, U(a, b, z) is the Tricomi's confluent hy- 
pergeometric function well-behaved at infinity, z — > oo, 
for negative E, and a = 1/2 — (E =p pBB)/haj c (=F cor- 
responds to spin up/down, respectively). Neglecting a 
small diamagnetic correction (quadratic in B) yields 



E = — €q + hw c m/2 ± /is-B, 



(3) 



where w c = eB/mt (mi, is the band mass), m = 
0,±1,±2, ... is the magnetic quantum number of the lo- 
calised state so that 
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a= \ + — 



(4) 



Hereafter we measure the magnetic field, b = B/Bq, in 
units of Bq = hn 2 /2e, where k = (2mb€o) 1 ^ 2 /h is the 
inverse localisation length of the zero-field state with the 
ionisation energy eo. 

Now using an integral representation of the confluent 
hypergeometric function [271 ] one obtains after integrat- 
ing over p in Eq. (TIJ 



II. HOPPING MAGNETOCONDUCTANCE VIA 
NON-ZERO ANGULAR MOMENTUM ORBITALS 

In organic and inorganic doped insulators lattice de- 
fects such as vacancies, interstitials, and excess neutral 
atoms or ions often localise carriers in finite momentum 
states rather than in the zero-momentum s-state [25j . 
Also the conventional nonmagnetic donors and accep- 
tors have nonzero orbital momentum states along with 
s-states, which are accessible for the hopping conduction 
at elevated temperature. Here we briefly outline the the- 
ory of hopping MR via non-zero angular momentum or- 
bits [in. 

In the hopping regime with localized carriers MR is 
caused by a strong magnetic field dependence of the expo- 
nential asymptotics of bound state wave functions at a re- 
mote distance from a donor (or an acceptor) . The asymp- 
totics is found using the Green function, G(r, r'; E), (GF) 
of free (or Bloch) electrons in the magnetic field with a 
negative binding energy E Q. The three-dimensional 
GF, G3_o(r, r'; E) is readily calculated using the Fourier 
transformation of the two-dimensional G 2 d (~$, p'\ E) 

m , 



G 3D (r,r';E) = JL J°° dpe^ z '-^G 2D {^,^';E- 
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where <f> and </>' are azimuth angles of and p ' respec- 
tively, I = (H/eB) 1 / 2 is the magnetic length, T(a) is the 
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Expanding the exponent in the square brackets in Eq. ([5]) 
up to the third power in x and performing the integration 
by the saddle-point method we obtain the asymptotics of 
the 3D impurity wave function, ip(r,B) cx G3£>(r,0; E) 
at 1 « nr < lib as [24| 



ip m (r,B) oc ip m (r,0) 
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(6) 

with r 2 = p 2 + z' 2 . 

For the s-wave bound state with m = this is the text- 
book asymptotics 0, Q accounting for the conventional 
positive MR quadratic in small B . On the contrary, for 
orbitals with nonzero orbital momentum the wave func- 
tion, Eq. (J6j , is linear in small B. Remarkably a local- 
ized state with a positive m expands in the magnetic field 
due to the magnetic lowering of its ionization energy by 
Huj c m/2, while states with negative ms shrink. The lin- 
ear term dominates in a wide range of realistic impurity 
densities for any nonzero m because of the small numer- 
ical factor (1/96) in the quadratic term in the exponent 
o'fEq.©. 

The hopping integral is proportional to ip m (r,B), 
where r is the distance between two hopping sites, which 
is assumed to be large compared with the localisation 
length. The hopping conductance is proportional to the 
hopping integral squared. If there is no time-reversal 
symmetry breaking, the states with the opposite direc- 
tion of the orbital angular momentum, m and — m, are 
degenerate, so that the linear term in the conductivity, 
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FIG. 1: Upper panel: Magnetoresistance to hopping via p- 
orbitals in TRS systems versus the magnetic field. Lower 
panel: The same MR in systems with the broken TRS. 



o-m cancels, and 



a{B) = a(0) 



nrb/2 



sinh(Kr&/2) 



cosh(m K r&)e- K ' Vrfc2 / 48 . 



(7) 

In this case the hopping conductivity, cr(B) first increases 
quadratically with the magnetic field (parabolic NMR) 
and only then decreases with B (positive MR), if up 2 /r < 
24m 2 — 4 as illustrated in FigfT] (upper panel) representing 
MR = a(0)/a(B) - 1 for p-states (to = ±1). 

Ions that carry a magnetic moment break the time- 
reversal symmetry and split to and —to states. Such 
zero-field splitting (ZFS) gives preference to hopping 
via orbitals with a lower ionisation energy (positive to). 
Hence in a ferromagnet with the frozen magnetisation 
MR for hopping via nonzero momentum orbitals is highly 
anisotropic changing from linear and negative in the field 
applied parallel to the magnetisation to linear but posi- 
tive in the opposite field. On the other hand if the global 
or local internal magnetisation rotates with the external 
magnetic field (superparamagnetism) magnetoresistance 
is negative and linear in small B no matter what the di- 
rection of the external field is. When the splitting due to 




FIG. 2: (Color online) Magnetoresistance measured at room 
temperature in ITO / PEDOT / AlqS/Ca device at bias volt- 
ages 12, 13, 14, 15, 17 V (symbols from bottom to top, re- 
spectively, Ref. \±M) described by Eq.fTUJ (solid lines) with 
Bh = 30 mT, m = 1 and the resistance ratio R shown in 
FigEl 



an internal (exchange) field is comparable or larger than 
the temperature, the hopping MR is found as 



MRh 
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sinh(Kr6/2) 
Krb/2 



exp 



K 3 p 2 r- 
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MRh in systems with the broken TRS, Eq.(j8|), shown in 
Fig Q] (lower panel), is quite distinguishable from MRh in 
TRS systems, FigQJ upper panel). 



III. OMAR 

Eq.([H]) allows us to address puzzling experimental ob- 
servations of negative NMR in a number of organic ma- 
terials [l8l - [2"(il | . There is experimental evidence for para- 
magnetic centers and ZFS in polymers, in particular 
in Alq3 j28j. More recently dilute magnetic impurities 
and magnetic domains have been observed in some tt- 
conjugated polymers [2{|. Hence as suggested in Ref. [13], 
if TRS is broken by such magnetic centers, the low-field 
magnetoresistance is dominated by the linear NMR via 
non-zero angular momentum orbitals, Eq.((5|). Here we 
extend our original description [24J of OMAR to the 
whole range of magnetic and electric fields used in the 
experiments. 

The device fabrication fl8l - [20l j started with glass sub- 
strates coated with a metallic layer, indium-tin-oxide 
ITO or the conducting polymer PEDOT. The semicon- 
ducting polymer (e. g. Alq3) layer was thermally evap- 
orated onto the bottom electrode, yielding an organic- 
semiconductor layer thickness of about 100 nm. The 
cathode, either Ca with an Al capping layer, Al, or Au 
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B (mT) 



was then deposited on top. 

The resistance of the device is the sum of the hop- 
ping resistance, Rh, of the organic semiconductor layer 
probably including the interface, and of other layers and 
contacts, R m . Hence, the magnetoresistance is expressed 
via metallic/contact MRm = R m (B)/R m (0) — 1 and the 
hopping magnetoresistance, MRh of Alq3 as 



MR: 



MRm MRh 



1 + 1/R 1 + R' 



(9) 



where R = i? m (0)/i?^(0) is the ratio of the zero- field 
resistance of other layers and contacts to the zero-field 
organic-semiconductor resistance. Reducing the number 
of fitting parameters we further assume that MRm is 
small, MRm <C MRh, and the first term in Eq.@ can 
be neglected. Also as outlined above the quadratic term 
in the exponent of Eq.® can be dropped in the relevant 
magnetic fields which yields 



MR 



(2B h smh(B/2B h )/B) 2 exp(-\m\B/B h ) - 1 
1 + R 



(10) 

where Bh = Bo/nr is the characteristic scaling field. 

As shown in Figl2l Eq.(fT0|) with R as a single fitting 
parameter describes remarkably well the experimental 
NMR for all magnetic and electric fields used in the ex- 
periment. The resistance ratio R increases with the volt- 
age, V, on the device, FigEl presumably due to a drop 
of the zero-magnetic field Rh(0) with V. 



IV. GRAPHITE MR 

Ref.[23[ observed unusual MR in commercially avail- 
able HOPG samples with different mosaicity and the 
room temperature out-of-plane/basal plane zero-field re- 
sistivity ratio pel Pab as high as 10 5 . The magnetic field 
was applied parallel to the hexagonal c-axis {B || c), and 
the in- plane p a b{B, T) was measured placing silver pasted 
electrodes on the sample surface. All resistance measure- 
ments were performed in the Ohmic regime (more sample 



FIG. 4: (Color online) Giant in-plane magnetoresistance of 
HOPG in the metallic (Y) direction at different temperatures 
(symbols [H) fitted with (B/B ln ) A/3 (lines). 



details and complementary magnetization measurements 
are found in Ref . p3l| ) . 

The graphite samples described in Ref. [23[ show quali- 
tatively different temperature dependence of the in-plane 
zero-field resistance for the current in two perpendicu- 
lar in-plane directions. In one direction (called here Y) 
the resistance is metallic-like while in the other direction 
(X) it is insulating-like in a wide temperature interval 
from 2K to the room temperature. Since a sufficiently 
strong magnetic field of about 1 Tesla or more makes 
graphene planes electrically isotropic the anisotropy is 
of electronic origin. The in-plane electrical anisotropy is 
most probably associated with an inhomogeneous carrier- 
density distribution, such that well doped metallic clus- 
ters are partially overlapped in Y-direction while they 
are separated by poorly doped insulating regions in X- 
direction. The in-plane anisotropy has been observed 
only in most anisotropic (p c / Pb ~ 10 5 ) graphite samples 
pointing to its quasi-2D origin. Our measured Kish and 
natural graphite crystals possess much lower resistivity 
ratio (p c / Pb ~ 10 2 -10 3 ) with no such effect observed. No 
in-plane anisotropy has been observed in HOPG samples 
with pd Pb < 10 4 either. It brings us to a conclusion that 
the electrical anisotropy is closely related to the reduced 
electron dimensionality in quasi-2D graphite providing a 
quasi-lD percolation in the planes. 

When a relatively weak magnetic field is applied per- 
pendicular to the planes, the MR along the metallic Y 
direction appears huge and positive, FigfJJ A giant pos- 
itive MR is naturally expected in doped graphite with 
virtually massless Dirac fermions [3(| since the param- 
eter {3 — uj c t becomes large already in the mT-region 
of the field (here r is the scattering time). However, 
it is neither quadratic as in the Boltzmann theory, nor 
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FIG. 5: Characteristic field B in of metallic clusters in HOPG 
as a function of temperature. 



linear in B as in the quantum magnetoresistance [Io| . 
It has been suggested [23[ that inhomogeneities are re- 
sponsible for the strong departure from these regimes. 
They lead to a radical rearrangement of the current flow 
pattern changing the magnetic field dependence of the 
transverse conductivity. Importantly, when /3 > 1 , then 
even relatively small inhomogeneities in the carrier den- 
sity lead to the MR proportional to B 4 ! 3 [3l|. In fact, 
MRy = (B / Bi n ) 4 / 3 perfectly fits the observed magnetic 
field dependence of MR in the metallic current direction, 
Fig.® with a single scaling parameter Bi n depending 
on fluctuations in the carrier density [3l| and tempera- 
ture, Fig(5] The temperature dependence of B in (T) is 
rcminicent of the temperature dependence of metallic re- 
sistivity, as it should be since B in oc 1/r 31]. 

The MR in the insulating direction is also anomalous, 
FiglU It is linear and negative in very low fields and pos- 
itive and superlinear in higher fields above the crossover 
point. The insulating- like temperature dependence of the 
zero-field resistance [23[ in this direction supports the 
view that the metallic clusters are virtually nonoverlaped 
along X, and the resistance is the sum of the hopping re- 
sistance, Ru of the insulating regions with low doping, 
and the metallic- like resistance R rn . Hence, the magne- 
toresistance in X direction is expressed via metallic MRy 
and the hopping magnetoresistance, MRh of insulating 
layers as in Eq.©, 



MRx 



MRy 



MRh 



1 + 1/R 1 + R' 



(11) 



where R = i? TO (0)/i?/i(0) is the ratio of the zero- field 
metallic resistance to the zero-field insulating resistance. 
This ratio ranges from about 0.1 at 2K to 1 at room 
temperature [23 1. 

However, different from OMAR Eq.(j9|) the first term 
in Eg. ([lip is signifincat and responsible for the upturn 
of MRx from negative to positive at some magnetic field 
about 70 mT, Fig.© as described in Ref . [23| . Here we 
analyse in more detail the low-field asymptotics of MRx. 
If TRS is broken in HOPG then according to Eq.fnj]) 
this asymptotics is linear in B because the first term in 
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FIG. 6: (Color online) Upper panel: Magnetoresistance in 
the insulating (X) direction (symbols). Lower panel: Low- 
field negative linear MR fitted as MRx = —B/B s with 
B s = (1 + R)Bh/\m\, shown in Fig[7] at different tempera- 
tures, indicating the time-reversal symmetry breaking. 



Eq. ([TTj) is superlinear, 



MRx : 



\m\B 



(1 + R)B h 



(12) 



for B -> 0. 

With increasing temperature the resistance ratio R in- 
creases and deeper localized states with a higher ioniza- 
tion energy become accessible for the hopping conduc- 
tance, so that the slope 1/Bs of the linear negative MR 
drops, as observed, Figs. ([6l7p . 



V. CONCLUSIONS 

In conclusion, we described the unusual magnetore- 
sistance of organic insulators (OMAR) [181420T ] and the 
highly- anisotropic in-plane magnetoresistance of quasi- 
two dimensional HOPG graphite [23J using the same the- 
ory of hopping magneto-conductance via non-zero angu- 
lar momentum orbitals [13]. While the theory with one 
or two scaling parameters provides accurate agreement 
with OMAR and the HOPG data, there are other the- 
oretical mechanisms of NMR outlined in the introduc- 
tion and beyond. Most of them explain NMR as some 
spin-correlation effects. Importantly, there is no NMR 
and virtually no MR in the same range of the magnetic 
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FIG. 7: Inverse slope B s of the linear NMR in HOPG versus 
temperature. 



field parallel to the planes, which rules out a spin ori- 
gin of the observed NMR in our quasi-2D HOPG sam- 
ples. Also there are no quantum magnetic oscillations 
at high temperatures, where unusual MRs are still ob- 
served, Figs. (l4l6[ ). so they are not related to the Landau 
quantization. In contrast to a number of spin correlation 
and weak localisation scenarios we observed the perfectly 
linear NMR at very low B in a wide temperature range, 
Fig.©, which according to our theory is a clear signature 
of the time-reversal symmetry breaking. 

The origin of the broken TRS in organic compounds 
and in graphite could be basically different. Instead 
of the (super)paramagnetism in polymers [HI, l29j the 
graphite samples of Ref. ^23] show a large diamagnetism. 
Naturally, some local paramagnetic centers responsible 



for the broken TRS in organic semiconductors could be 
also found in graphite, with their magnetic response over- 
whelmed by the large diamagnetism of metallic clus- 
ters. However the observed large diamagnetism could 
suggest a more intriguing mechanism of TRS breaking, 
such as superconducting cluste rs I3 2 | with an unconven- 
tional (chiral) order parameter [33l.l34|. There is a kink in 
the field dependence of the diamagnetic magnetization of 
the HOPG samples at B k ss 0.2 T, H3|, resembling the 
behavior of type-II superconductors in magnetic fields 
exceeding the lower critical field. Supporting this possi- 
bility the electrically inhomogeneous samples are becom- 
ing homogeneous insulators in sufficiently high magnetic 
fields, which could suppress the superconductivity. Also 
in some 2D lattices there is the broken time-reversal sym- 
metry (i.e. some internal magnetic ordering) in the nor- 
mal state with a spontaneous quantum Hall effect but 
without any net magnetic flux at high temperatures [35[ . 
At lower temperatures a pair of spontaneously generated 
current loops in adjacent graphene layers, having odd- 
parity with respect to the two layers, could also break 
TRS HU 

More generally our findings point to an inhomogeneous 
doping and a semiconducting gap in graphite. 
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